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Popular Summary 

Since first light in early 2000, operational global quantitative retrievals of aerosol 
properties over land have been made from MODIS observed spectral reflectance. These 
aerosol products have been used for dozens of applications, including some not even 
imagined before launch. We have been continuously evaluating and validating the 
algorithm and products, and have noted opportunities for their improvement. In early 
2006, we replaced the original algorithm. This new algorithm (known as Version 5.2 or 
V5.2) performs a simultaneous inversion of two visible channels (0.47 and 0.66 pm), plus 
a shortwave-IR (2,12 pm) channel, thus improving sensitivity to the scattering of coarse 
aerosols (> 1 pm in radius) such as dust. Also, we have improved the assumptions related 
to surface reflectance, aerosol optical properties, and radiative transfer calculations. 

Inversion of the three channels yields three nearly independent parameters, the aerosol 
optical depth (x) at 0.55 pm, the non-dust or fine weighting (q) and the surface 
reflectance at 2.12 pm. Retrievals of small magnitude negative x values (down to -0.05) 
are now considered valid, thus normalizing the statistics of x in near zero x conditions. 

On a “test-bed” of 6300 granules from both Terra and Aqua satellite platforms, the 
products from V5.2 show marked improvement over those from the previous versions. 

Compared with surface based sunphotometer (AERONET) observations, regression of 
MODIS x (at 0.55 pm) has an equation of: y = l.Olx + 0.03, R = 0.90. Mean x for the 
test-bed is reduced from 0.28 to 0.21, which is more in line with estimates from transport 
models. We also added some new products for V5.2, including “fine aerosol x” and 
spectral surface reflectance. The new algorithm began operational production of 
“Collection 005” in April 2006, and is expected to finish re-processing the entire MODIS 
mission by the end of 2006. 
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A new algorithm for retrieving aerosol properties over land from MODIS spectral 

reflectance 

Robert C. Levy, Lorraine A. Remer, Shana Mattoo, Eric Vermote, Yoram J. 

Kaufman 

Since first light in early 2000, operational global quantitative retrievals of aerosol 
properties over land have been made from MODIS observed spectral reflectance. These 
products have been continuously evaluated and validated, and opportunities for 
improvements have been noted. We have replaced the original algorithm by improving 
surface reflectance assumptions, the aerosol model optical properties and the radiative 
transfer code used to create the lookup tables. The new algorithm (known as Version 5.2 
or V5.2) performs a simultaneous inversion of two visible (0.47 and 0.66 pm) and one 
shortwave-IR (2.12 pm) channel, making use of the coarse aerosol information content 
contained in the 2.12 pm channel. Inversion of the three channels yields three nearly 
independent parameters, the aerosol optical depth (t) at 0.55 pm, the non-dust or fine 
weighting (rj) and the surface reflectance at 2.12 pm. Finally, retrievals of small 
magnitude negative x values (down to -0.05) are considered valid, thus normalizing the 
statistics of x in near zero x conditions. On a ‘test bed’ of 6300 granules from Terra and 
Aqua, the products from V5.2 show marked improvement over those from the previous 
versions, including much improved retrievals of x, where the MODIS/AERONET x (at 
0.55 pm) regression has an equation of: y = l.Olx + 0.03, R = 0.90. Mean x for the test 
bed is reduced from 0.28 to 0.21. 
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1. Introduction 


Aerosols are major players in Earth’s climate, radiation budget, cloud processes 
and air quality, and increasingly sophisticated and accurate remote sensing techniques 
have been introduced to understand aerosols and their effects. Especially for aerosols 
over land, the first operational global satellite dataset has been provided by the Moderate 
Imaging Spectrometer (MODIS). Since MODIS’ launch aboard Terra (in late 1999) and 
aboard Aqua (in early 2002), the use of the MODIS aerosol products has grown 
exponentially. Since launch, MODIS data and specifically aerosol data have been used 
for dozens of applications and in hundreds of publications. Not only have MODIS aerosol 
products been used to answer (intended) scientific questions about radiation and climate 
(e.g. IPCC, 2001; Yu et al., 2006), they are being used for applications not previously 
intended. One example is using MODIS to monitor surface air quality for health (e.g. 

Chu et al., 2003; Al-Saadi et al,, 2005). 

The operational algorithm over land has been using MODIS reflectance data in 
three channels (0.47, 0.66 and 2.12 pm; channels 3, 1 and 7) to retrieve total spectral 
‘aerosol optical depth’ (AOD or x) and ‘Fine aerosol Weighting’ (FW or r]), reported at 
0.55 pim. Since launch, the aerosol products have been monitored for quality, and the 
algorithm has been continuously updated. Details of a previous version (V4.2) and the 
products created for ‘Collection 004’ (C004) has been described in Remer et al., (2005). 
The last update to the algorithm was known as ‘V5.1\ but it never became operational. 

In order to be applied in both climate and air pollution applications, MODIS 
aerosol retrievals must meet certain expected accuracy (Kaufman et al. 1997a). MODIS 
should be able to retrieve x to within expected errors, specifically: 

Ar = ±0.05 ± 0. 15r (1) 

(Remer et al., 2005). To this end, a number of papers have attempted to ‘validate’ the 
retrieved properties of C004 and before, by comparing MODIS derived values to 
standard (ground truth) aerosol measurements, using the co-location theory of Ichoku et 
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al., (2002). Ground based sunphotometers, especially from the Aerosol Robotic 
NETwork (AERONET - (Holben et al., 1998)), have provided the bulk of the 
comparison data (e.g. Chu et al., 2002; Levy et al., 2005; Remer et al., 2005; Ichoku et al. 
2005). Most of these validation studies have shown that although MODIS generally 
derives x within the expected error, MODIS tends to over-estimate x for small x and 
underestimate for high x (Chu et al, 2002; Remer et al. 2005; Levy et al. 2005). In fact, at 
0.55 pm, the consensus is a MODIS/AERONET x regression of approximately 

T MODIS = 0.1 + 0.9 "C AERONET- (2) 

Why is the regression not one to one? The algorithm made a number of 
assumptions as to the aerosol optical properties, surface reflectance, MODIS channel 
central wavelengths, Rayleigh (molecular) optical depth, radiative transfer, transparency 
of 2.12 pm channel, and retrieval philosophy /logic. This paper introduces a new aerosol 
retrieval algorithm, known as ‘ V5.2’ that analyzes the validity of each assumption in the 
C004 algorithm. Section 2 introduces the C004 MODIS products and AERONET data 
used in this study. Section 3 summarizes the new aerosol models and look-up table 
described by Levy et al., (2006). The surface reflectance properties are discussed in 
Section 4. Section 5 introduces a new retrieval philosophy, and section 6 discusses the 
products that will be Collection 5 (C005). Finally, we show some provisional validation 
of V5.2 in Section 7. 

2. The MODIS C-004 and AERONET L2A datasets 

For this work, we made extensive use of both MODIS and AERONET data 
products. Aerosol products have been derived from Terra reflectance observations since 
2000 and since 2002 from Aqua. As of early 2005, most MODIS observations (‘Level 1’) 
through 2004 had been processed or re-processed into ‘Level 2 products’ (L2) using 
consistent retrieval algorithms. This collection of products is known as ‘Collection 4’ or 
‘C004’. At some sites, AERONET has been reporting since 1993, and as of early 2005, 
most of the AERONET data have been re-processed and quality assured by the 
AERONET team. These products are also known as ‘Level 2’, but to differentiate them 
from the MODIS products, we denote the AERONET products as ‘L2A’. 
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The MODIS instruments aboard Terra and Aqua both measure spectral radiance 
in 36 channels, in resolutions between 250 m and 1 km (at nadir). In polar-orbit about 
700 km above the earth, MODIS views a swath about 2300 km, resulting in near daily 
global coverage of Earth’s land/ ocean/atmosphere system. The swath is broken into five- 
minute ‘granules’, each about 2030 km long. The aerosol algorithm over land makes use 
of gas-absorption corrected spectral solar reflectance measurements (in visible through 
the IR bands) to perform cloud masking and pixel selection, and then retrieve aerosol 
optical depth (AOD or r/and fine-dominated aerosol fraction (known as ‘fine weighting’, 
FW or rj) at 10 km resolution (at nadir). 

The basic concepts of the land algorithm was introduced by Kaufman et al., 
(1997a) and updated by Remer et al,, (2005) to describe the products of C004. From pre- 
launch through C004, the algorithm has gone through a series of updates, detailed at 
http://modis-atmos.gsfc.nasa.gov/MQD04 L2/historv . The theory of the MODIS over- 
land algorithm is as follows. The upward spectral ‘reflectance’ (normalized solar 
radiance) at the top of the atmosphere (TOA) is a function of successive orders of 
radiation interactions within the coupled surface-atmosphere system. The TOA angular 
(6 0 ,6,<p = solar zenith, view zenith and relative azimuth angles) spectral reflectance 
( p x (0 o ,6,<p)) at a wavelength X results from: scattering of radiation within the atmosphere 
without interaction with the surface (known as the ‘atmospheric path reflectance’), the 
reflection of radiation off the surface that is directly transmitted to the TOA (the ‘surface 
function’), and the reflection of radiation from outside the sensor’s field of view (the 
‘environment function’). The environment function is neglected so that to a good 
approximation: 


i-s ;J o;.(0 o ,0,0) 


(3) 


where FdK is the ‘normalized downward flux’ for zero surface reflectance, T x represents 
‘upward total transmission’ into the satellite field of view, s x is the ‘atmospheric 
backscattering ratio’, and ff x is the angular ‘surface reflectance’. Except for the surface 
reflectance, each term on the right hand side of Equation 3 is a function of the aerosol 
type and loading (t). Assuming that a small set of aerosol types and loadings can describe 
the range of global aerosol, we have created a lookup table that contains pre-computed 
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simulations of these aerosol conditions. The goal of the algorithm is to use the lookup 
table to determine the conditions that best mimic the MODIS-observed spectral 
reflectance ff x , and retrieve the associated aerosol properties (including x and rj). The 
difficulty lies in making the most appropriate assumptions about both the surface and 
atmospheric contributions. 

The sunphotometers of the Aerosol Robotic NET work (AERONET) provide a 
comprehensive data set of aerosol properties. These include direct ‘sun’ measurements of 
spectral x in four or more wavelengths (to include 0.44, 0.67, 0.87 and 1.02 pm), and 
indirect ‘sky’ measurements that lead to estimates of aerosol optical properties and 
aerosol size distributions (Holben et al., 1998). These data go through rigorous 
calibration and cloud screening processes, resulting in the L2A products. The AERONET 
direct sun measurements are made approximately every 15 minutes during mid-day and 
more often during sunrise and sunset. In addition to spectral x, AERONET also provides 
estimates of columnar water vapor w (precipitable water in units of [cm]). O’Neill et al. 
(2003) developed a method for estimating r] from the direct sun measurements of spectral 
x. The AERONET sky radiance measurements are made less often (about once per hour), 
and are inverted simultaneously either assuming spherical aerosol particles (Dubovik and 
King, 2000) and/or spheroid particles (Dubovik et al., 2002; Dubovik et al., 2006). Under 
either particle assumption, the fundamental derived parameters include spectral x, 
spectral complex refractive index, the volume distribution as a function of 22 radius size 
bins (dV/dlnR), and fitting error to the radiance measurements. Additional parameters are 
then calculated that include Angstrom exponents, properties of two (fine and coarse 
mode) lognormal aerosol distributions, spectral single scattering albedo (SSA or a > 0 ) and 
asymmetry parameter (g) of the lognormal modes. 

Although the actual products provided by MODIS and AERONET are not 
necessarily physically identical, in many cases they are comparable. For example, 

MODIS retrieves x at 0.55 pm, whereas AERONET retrieves at 0.44 and 0.67 pm, and 
some instruments also measure at 0.50 pm. However, fitting a quadratic equation through 
the logarithms of both x and wavelength, AERONET x can be interpolated to 0.55 pm 
(Eck et al., 1999). Comparison of q is trickier. Over land, MODIS considers q to be the 
contribution of the fine-dominated model (the non-dust model) to the total x, the 
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AERONET sky retrievals designate rj to be the volume contribution from aerosol below a 
radius of 0.6 pm, whereas the O’Neill method separates fine and coarse aerosol by 
spectral behavior. Practically, however, the definitions of r] are similar enough so that 
they should be correlated (Kleidman et al., 2005, Anderson et al., 2005, Chu etal. 2005). 

Over 15,000 pairs of MODIS and AERONET ‘ sun’ data, at over 200 global sites, 
have been co-located in time via the technique of Ichoku et al., (2002). A valid 
MODIS/AERONET match is considered when there at least five (out of a possible 25) 
MODIS retrievals (10 km x 10 km resolution) within the box, and at least two (out of a 
possible five) AERONET observations within the hour. This co-located data set was used 
for a number of applications described in this document, including studies of surface 
reflectance, and validation of MODIS products. 

3. New aerosol lookup tables 

A number of studies (e.g. Chu et al., 2002, Reiner et al., 2005, Levy et al., 2005) 
have demonstrated that MODIS/AERONET regression of t over land results in slope less 
than one; meaning that MODIS tends to under-estimate optical depth, especially as the 
optical depth increases. Ichoku etal., (2003) and Levy et al., (2005) found that updating 
the assumed aerosol properties in Southern Africa and the U.S. East Coast, respectively, 
improved the retrieval in those regions. Results from studies such as these led Levy et al., 
(2006) to consider deriving new aerosol optical models for V5.2, using over 136,000 
AERONET sky retrievals. 

Levy et al., (2006) performed a cluster analysis of spherical AERONET 
almucantur inversions, and found that global fine-mode dominated aerosol could be 
separated into three types, varying primarily by single scattering albedo (SSA or co 0 ). 
There was a ‘non-absorbing’ aerosol model (co 0 ~0.95) that presumably corresponded to 
urban aerosol in the industrialized northern hemisphere, an ‘absorbing’ (co o ~0.85) aerosol 
model found in the known savanna-burning regions of South America and Africa, and a 
‘neutral’ aerosol model found in primarily forest fire regions and the developing world. 

A similar analysis of spheroid retrievals showed that only a single model was necessary 
to describe the properties of presumably dust aerosol. They then ‘fixed’ the aerosol type 
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at a given season and location, based on the dominant aerosol type found during 
clustering. These decisions were mapped onto a 1° longitude x 1° latitude grid, such that 
a fine aerosol type is assumed for each grid point, globally. This global map approach, 
that is not hardwired into the processing code, allows for easy alterations as new 
information becomes available. 

For the set of new assumed spherical fine aerosol models (absorbing, non-absorbing 
and neutral), Levy et al., (2006) computed the V5.2 LUT using a combination of MIEV 
(Wiscombe et al., 1980) and RTS (Evans and Stephens, 1991). For the non-spherical 
(assumed spheroids) coarse (dust) aerosol they used the T-matrix code described in 
Dubovik et al., (2002, 2006). In addition, corrections for Rayleigh optical depths and 
center wavelengths were performed. The LUT was calculated for seven aerosol loadings 
( r o .55 = 0-0, 0.25, 0.5, 1.0, 2.0 3.0, and 5.0), 9 solar zenith angles (6 0 = 0.0, 6.0, 12.0, 24.0, 
35.2, 48.0, 54.0, 60.0 and 66.0), 16 sensor zenith angles (0 = 0.0 to 66.0, increments of 
6.0), and 16 relative azimuth angles ( <p = 0.0 to 180.0 increments of 12.0). For specific 
AERONET sites, sun-derived spectral x dependence were compared with the assumed 
models. 

4. VIS/S WIR surface reflectance assumptions 

When performing atmospheric retrievals from MODIS or any other satellite, the 
major challenge is separating the total observed reflectance into atmospheric and surface 
contributions (e.g. Equation 3), and then defining the aerosol contribution. Over the 
ocean, the surface is nearly black at red wavelengths and longer, so that assuming 
negligible surface reflectance in these channels is a good approximation. Over land, 
however, the surface reflectance in the visible and SWIR is far from zero and varies over 
surface type. As the land Surface and the atmospheric signals are comparable, errors of 
0.01 in assumed surface reflectance will lead to errors on the order of 0.1 in x retrieval. 
Errors in multiple wavelengths can lead to poor retrievals of spectral x, which in turn 
would be useless for estimating size parameters. 

Kaufman and colleagues (e.g. Kaufman et al., 1997b) observed that over vegetated 
and dark soiled surfaces, the surface reflectance in some visible wavelengths correlated 
with the surface reflectance in the SWIR. Parallel simulations by vegetation canopy 
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models, showed that the physical reason for the correlation was the combination of 
absorption of visible light by chlorophyll and infrared radiation by liquid water in healthy 
vegetation (Kaufman et al., 2002). These relationships were such that the surface 
reflectance values in the visible were nearly fixed ratios of that in the SWIR (Kaufman et 
al., 1997b). As applied for C004, surface reflectance at 0.47/rm (channel 3) and 0.66 /<m 
(channel 1) were assumed to be one-quarter and one-half, respectively, of the surface 
reflectance in the mid-SWIR 2.12/rm (channel 7) (Kaufman et al., 1997b). 

Regression of C004 (and prior) MODIS-derived x to AERONET sunphotometer data 
(Chu et al., 2002; Remer et al., 2005) has shown that while the products generally agreed 
(~60-65%) to within the expected errors of Eq. 1 (±0.05 ± 0.15x), there was a positive 
offset of about 0.1 (Eq. 2). This means that the C004 algorithm generally over-estimates x 
in pristine conditions, which questions the assumed surface reflectance, From data 
observed during the CLAMS experiment of 2001, Levy et al., (2005) found that higher 
values of VIS/SWIR surface ratios (e.g. 0.33 and 0.65 for the blue and red, respectively) 
improved the continuity of the MODIS over-land and over-ocean aerosol products along 
the coastline of the DelMarVa Peninsula. The MODIS/AERONET x regression over 
near-coastal sites was also improved. However, at locations far from the coastline, the 
CLAMS VIS/SWIR ratios tended toward over-correction of the surface reflectance and 
retrievals of x less than zero. Thus, we know that a single ratio is not globally applicable. 

It is known that earth’s surface is not Lambertian, and that some surface types exhibit 
strong bi-directional reflectance functions (BRDF). Gatebe et al., (2001) flew the Cloud 
Absorption Radiometer at low altitudes over vegetated surfaces and found that the 
VIS/SWIR surface ratios varied as a function of angle, and often greatly differed from the 
one-quarter and one-half ratios assumed in C004. Remer et al., (2001) also noted that the 
VIS/SWIR surface ratios varied as a function of scattering geometry. In fact, under 
certain geometry, the assumed VIS/SWIR surface relationship broke down completely. 

To continue with this philosophy of using VIS/S WIR relationships to determine VIS 
surface reflectance, we will need to link VIS/SWIR relationships to geometry and to 
surface characteristics that we can measure from space without confusion from overlying 
aerosol layers. Note that alternative philosophies including using digital surface type 
models and global maps of measured spectral albedo were explored, but found to be less 
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useful in the aerosol retrieval than measured parameters and empirical relationships. 
However, in order to develop the empirical relationships we need a data base of surface 
reflectances that are representative of a wide range of global conditions. Before Terra 
launch such data were unavailable. Today, we can use the accumulation of MODIS data 
to derive a large database of surface reflectances. 

4.1. Atmospheric correction of C-004 MODIS/AERONET co- 
located products 

Atmospheric correction (Kaufman and Sendra, 1988) attempts to calculate the 
optical properties of the surface, by theoretically subtracting the effects of the atmosphere 
from the satellite-observed radiation field. One needs to assume the optical properties of 
the intervening atmosphere, including all aerosol and non-aerosol components. In 
addition to knowing or assuming all atmospheric components, accurate radiative transfer 
(RT) is also required. The atmospherically corrected surface reflectance is calculated 
by re-arranging Eq. 3. 

In order to minimize errors arising from multiple scattering by the aerosol, we 
have limited our atmospheric corrections to conditions of x in the green less than 0.2. Out 
of the original 15,000 co-located MODIS/AERONET points (described in section 2), 
there are over 10,000 collocations with low x. The archive includes the ‘gas absorption 
corrected’ MODIS-Level 2 observed reflectance, as well as AERONET-observed (L2A) 
spectral x and column water vapor depth. For atmospheric correction we use the same 
reflectances that the MODIS algorithm used to retrieve the x in the 10 km box that 
contains the AERONET instrument. These reflectances have been corrected for gas 
absorption, and clouds, snow, inland water and bright surfaces have been eliminated. , 
Aerosol and water vapor characteristics for the correction are provided by AERONET 
observations in that 10 km box. The AERONET data is average over an hour, centered at 
over pass time. The molecular properties of the atmosphere are assumed those of the U.S. 
standard atmosphere, with the Rayleigh optical depth (ROD) values scaled from sea level 
values, -according to the elevation/air pressure of the sunphotometer. 
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The relation between the satellite-measured reflectance and the surface 
reflectance is a complicated function of the atmospheric effects of scattering and 
absorption by the aerosol. Previous atmospheric correction exercises often assumed some 
form of the Continental aerosol model (e.g. Vermote et al. 1997), to describe both the 
scattering and absorption properties. While this model may provide reasonable 
simulations in blue and red wavelengths, it cannot be expected to provide accurate 
simulations in the MODIS 2.12 pim channel, especially for either extreme of fine or 
coarse -dominated aerosol conditions. Instead, we used AERONET-derived Angstrom 
exponent (a) to decide whether the scene was fine or coarse aerosol. The ‘neutral’ 
(generic/developing world; SSA ~ 0.9) aerosol type (Levy et al. 2006) was assumed for 
the 4200 cases where col. 6. Values of a< 0.6 led to assuming the coarse (dust) model 
(about 400 cases). Co-locations where 0.6<a<1.6 (about 6000 cases) were not used due 
to uncertainties of aerosol mixing. 

4.2. Mean values of VIS/SWIR surface reflectance relationships 

Atmospheric correction was performed on the 4600 MODIS/AERONET co- 
locations having r 0 . 55 < 0.2 and a<0.6 or col.6. Figure 1 plots the entire set of 
atmospherically corrected visible surface reflectance (in the blue ff 0J7 and the red ff 066 ) 
versus that in the mid-SWIR and their regression lines. While not plotted, also 
considered were the regressions if they were forced through zero, thereby assuming that 
zero SWIR reflectance is zero reflectance over the entire spectrum (which would be 
equivalent to deriving simple ratios). Correlation (R) values are 0.93 for the red, but only 
about 0.75 for the blue. In the blue, forcing a regression through zero is quite different 
than that not constrained. If forced through zero, the slope tends toward 0.36, whereas 
including the offset (about +0.01 1) yields a slope closer to the assumed one-quarter 
(0.258). In the red, whether including offset or not, the slope is about 0.55. Thus in a 
mean sense, atmospheric correction of MODIS data yields VIS/SWIR surface reflectance 
relationships that differ substantially from the assumed C004 VIS/SWIR ratios. Fitting 
blue to red (Figure 8b) has higher correlation and less scatter than blue to SWIR, 
specifically R = 0.87. There is less difference between fitting through zero and not, such 
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that a straight blue/red ratio is about 0.54, and the full regression has slope = 0.508 and 
offset = 0.008. Therefore, instead of the 0.47 pm and 0.66 |xm surface reflectance being 
calculated separately from 2.12 pm, we will calculate the 0.66 pm surface reflectance 
from that in 2. 12 pm, followed by calculating 0.47 pm from the 0.66 pm, i.e. 

Po .66 ~ /(P2.12) ^ 

P0.47 = S(P 0.66^ 

where f() and g() are two independent relationships. To test the robustness of the 
relationship we performed similar regression with only the 2058 points where r O J5 <0.1. 
The relationship stayed nearly the same, except with slightly higher correlation values. 

4.3. Variability of VIS/SWIR surface reflectance relationships: 
Angle 

As noted in Fig. 1 the VIS/SWIR surface reflectance regressions display large 
scatter. For example, where 2.12 pm surface reflectance is 0.15, simply assuming the 
mean values of the red/SWIR and blue/red relationships would result in estimating the 
0.66 pm surface reflectance as 0.083 and 0.47 pm surface reflectance of 0.050. The 
scatter plots show that in reality, the 0.66 pm reflectance could vary between 
approximately 0.05 and 0.1, and the 0.47 pm surface reflectance between 0.01 and 0.07. 
Obviously, this could result in very large errors in retrieved x, on the order of 0. 1 or 0.2 
or more. Therefore, to reduce the scatter we look for dependencies on other parameters to 
refine the relationships. 

Gatebe et al. (2002) and Remer et al. (2001) suggest that the VIS/SWIR surface 
reflectance relationships are angle dependent. Out of different possible angle parameters 
(solar zenith angle, sensor zenith angle, glint angle or scattering angle) we found that the 
scattering angle had the largest influence on the VIS/SWIR surface reflectance 
relationship. The scattering angle, 0 , is defined as: 

0 = cos“ 1 (-cos0 o cos0 + sin6> o sin0cos^) (5) 
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where 6^9 and <p are the solar zenith, sensor view zenith and relative azimuth angles, 
respectively. The data from Fig.l were sorted according to scattering angle and put into 
20 groups of equal size (about 230 points for each scattering angle bin). Fig. 2 (a) 
displays the median values of surface reflectance in each bin as a function of scattering 
angle, and shows a definite relationship at 2.12 pm, less at 0.66 pm, and nearly none at 
0.47 p,m. Since Fig. 1 showed that both a slope and y-offset was necessary to regress VIS 
and SWIR surface reflectance, we look for scattering angle dependence on both 
parameters. Fig 2 (b-d) plots the slope, y-offset and correlation of the surface reflectance 
relationships calculated in each scattering angle bin and plotted as a function of scattering 
angle. The P S 0.66 / P S 2.12 regression slope (r0660 in the figure) shows dependence on 
scattering angle, whereas the p s 0 . 47 /p s 0 .66 regression slope (rVIS in the figure) shows 
nearly none. The regressed y -intercept shows strong dependence on scattering angle for 
both relationships. Especially interesting is that the red/SWIR y-offset goes from 
positive to negative with increasing scattering angle, with a value of zero near 0=135°. 

4.4. Variability of VIS/SWIR surface reflectance relationships: 
Surface type and MV1 

Because AERONET sites are located in different surface type regimes, it could be 
expected that the VIS/SWIR surface relationships will vary based on surface type and/or 
season. Using the International Geosphere/Biosphere Programme’s (IGBP) scene map of 
USGS surface types and formatted for MODIS validation 

(http://edcdaac.usgs.gov/modis/modl2clv4.asp), we determined the scene type of the 
MODIS/AERONET validation box. We then separated urban from non-urban surfaces, 
and grouped into season (winter or summer) and general location (mid-latitude or 
tropical). Generally, “greener” surfaces (midlatitude summer sites both urban and 
nonurban) have higher red to SWIR ratios (red/ S WIR>0. 55) than winter sites 
(red/S WIR<0.55) or tropical savannas and grasslands. As for the blue to red channel 
surface reflectance relationships, except for the urban sites during summer (blue/red ratio 
~ 0.766), the relationships around the globe are relatively consistent (blue/red ~ 0.52). 
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Except for urban areas, most surfaces seem to have VIS/SWIR surface reflectance 
relationships that vary as a function of their “greenness.” Can we relate the surface 
reflectance relationships to a vegetation index (VI)? The Normalized Difference 
Vegetation Index (NDVI), defined as a function of the red (0.66 pm - channel 1) and 
near-IR (0.86 pm - channel 2), can be heavily influenced by aerosol (Reference). An 
alternative is the NDVI swm , defined as: 

NDVI swm — (Pi 24 ~ P 2 .I 2 ) /(Pl .24 + P 2 .I 2 ) (6) 

where p 124 and p 2I2 are the MODIS-measured reflectances of the 1.24 pm channel 
(MODIS channel 5) and the 2. 12 pm channel (channel 7), which are much less 
influenced by aerosol (except for heavy aerosol or dusts). This index is also known as 
NDVI um (Mid-InfraRed) in the work of Miura et al. (1998) and others. In aerosol free 
conditions NDVI swm is highly correlated with regular NDVI. A value of NDVI swm > 0.6 is 
a highly vegetated area, whereas NDVI swm < 0.2 is representative of sparse vegetation. 
Figure 3 plots the relationship of the 0.66 pm channel and 2. 12 pm channel 
(atmospherically corrected) surface reflectance relationship, for nonurban sites, as a 
function of low, medium and high values of NDVI swm . Clearly, as the NDVI swm 
increases, the ratio between 0.66 pm and 2. 12 pm surface reflectance increases, and we 
will use this relationship in the final VIS/SWIR surface reflectance parameterizationA 

4.5. Final parameterization of VIS/SWIR surface reflectance 
relationships 

Results of the global atmospheric correction exercise imply that not only do the 
VIS/SWIR surface relationships differ from the ratios assumed by C004, they also have a 
strong dependence on both geometry and surface type. The new (V5.2) VIS/SWIR 
surface reflectance relationship is parameterized as a function of both NDVI swm and 
scattering angle 0, such that Equation 4 can be expanded: 
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Po.66 " f (P2.I2) — P2.12 * S l°P e 0.66 1 2.12 + .66/2.12 

and (7) 

Po.47 = 8 (Pom) ~ P0.66 * S l°P e OA 7 l 0.66 + yl n ^OA 7 IOj 66 

where 

slope 0( 66 / 2 . 12 = sloped™* + 0.0020 
3 ^ 0 . 66 / 2.12 = 0-000250 + 0.033, 
slope 047/066 = 0.49, and 
3^0.47/0.66 = 0.005 

where in turn 

slope N ™lll -0.48;A IDVI mR <0.25, 
slopeZIni = 0.58;A® VI SWIR > 0.75 (9) 

=0.48 + 0.2(NDVI swm - 0.25);0.25 £ NDVI swm s0.75 

Note that while the above parameterization was based on the results of Figs 1-3, the 
coefficients are not identical to those in the figures. The atmospheric corrected data set is 
the broadest and most comprehensive representation of global surface reflectance 
relationships, still it is limited to AERONET site locations, which are in turn are mostly 
concentrated in certain geographical regions. Trial and error was used to modify the basic 
results from the AERONET-based atmospheric correction, to give more realistic MODIS 
retrievals globally, (especially in places were few or.no AERONET sites are located). 

4.6. Notes on VIS/SWIR surface reflectance relationship errors 

We note that even with the surface reflectance parameterization, there still will be 
errors in estimating surface reflectance. According to the MODIS Land Surface 
Reflectance Homepage (http://modis-sr.ltdri.org/html/prodacc.html . improper aerosol 
model assumptions can lead to errors in atmospherically corrected reflectance on the 
order of 0.002 in the 0.47 and 0.66 pm channels, and 0,006 at 2.12 pm. The errors are 
especially large at 2.12 pm due to potentially choosing a fine-dominated model instead of 
a coarse-dominated model (or vice-versa). However, since our study pre-determined the 
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choice of fine or coarse-dominated aerosol models via the AERONET-observed 
Angstrom exponent, presumably errors at 2. 12 pm should be much less. The difference in 
spectral optical thickness between the V5.2 fine-dominated neutral and absorbing models 
(atx 055 =0.5) is about 0.02, 0.02 and 0.002, respectively in the 0.47, 0.66 and 2.12 pm 
channels. On average, this would be equivalent to errors of 0.002, 0.002 and 0.0002, 
respectively in surface reflectance, but would vary according to the differences in phase 
function. Regardless, the error at 2. 12 pm is small enough so that the derived surface 
reflectance relationship should be reasonably robust. 

Of course, other errors may creep into the surface reflectance parameterization. These 
include, but are not limited to additional surface BRDF effects lost during averaging over 
scattering angle and errors due to MODIS instrument calibration. The MODIS Land 
Surface Reflectance Homepage suggests that these errors can cause reflectance errors that 
similar in magnitude to those caused by improper aerosol model assumptions. 

5. Inversion of spectral reflectance: The V5.2 algorithm 

A major limitation of the old C004 algorithms is that aerosol is assumed 
transparent in the 2.12 pm SWIR channel. The surface reflectance in 2.12 pm is assumed 
to be exactly the value of the observed TOA reflectance in that channel. Under a dust 
aerosol regime, aerosol transparency is an extremely poor assumption. Even in a fine 
aerosol dominated regime, x is not zero. For example, for our “generic/developing world” 
(neutral co 0 ~0.90) aerosol, x 055 of 1.0 corresponds to x 212 of 0. 1 14. For a given angle (say 
0 O = 36°, 0 = 36°, and (J> = 72°) assuming x 212 = 0.0 instead leads to error in 2. 12 pm path 
reflectance of about 0.012. Via the VIS/SWIR reflectance relationship, the reflectance 
error at 0.66 pm would be on the order of 0.006, leading to ~ 0.06 error in retrieved x. As 
a percentage of the actual x, the error is not very large. However, combined with errors at 
0.47 pm, the resulting incorrect Angstrom exponent leads to error in estimating rj. 

In the spirit of the MODIS aerosol over ocean algorithm (Tan re et al., 1997), we 
developed a multi-channel reflectance inversion for retrieving aerosol properties over 
land. Analogous to the ocean algorithm’s combination of fine and coarse aerosol modes, 
the V5.2-land algorithm attempts to combine fine-dominated and coarse-dominated 


16 



aerosol models (each composed of multiple modes) to match with the observed spectral 
reflectance. The 2. 12 pm channel is assumed to contain both surface and aerosol 
information, and the visible surface reflectance is a function of the new Y5.2 VIS/SWIR 
surface reflectance relationships. Simultaneously inverting the aerosol and surface 
information in the three channels (0.47 pm, 0.66 pm and 2.12 pm) yields three 
parameters: r a55 , r/ 055 and the surface reflectance (/7 2J2 ). 

We rewrite equation 3, but note that the calculated spectral total reflectance p x at 
the top of the atmosphere is the weighted sum of the spectral reflectance from a 
combination of fine and coarse -dominated aerosol models, i.e. 

Pi = VP*/ + (1 - ri)Pi ( 10 ) 

where p* f x and p c x are each composites of surface reflectance ff x and atmospheric path 
reflectance of the separate aerosol models. That is: 

p'/ = P a /+F&T'pl/(l-s{pt) 

and ( 11 ) 

X-sIpD 

where p af x and p ac x are the fine and coarse model atmospheric path reflectance, F f dX and 
are normalized downward fluxes for zero surface reflectance, 1 f x and T x represent 
upward total transmission into the satellite field of view, and sf x and s c x are atmospheric 
backscattering ratios, The weighting parameter, r\ in Eq. 10, is defined for X = 0.55 pm. 
In the appendix of Remer et al. (2005) they show how this parameter also represents the 
fraction of the total optical thickness at 0.55 pm contributed by fine model aerosol. Note 
the angular and x dependence of some of the terms: p a =p a ( % 6 0 , 6,<p), F=F( r,6 0 ), 
T=T(r,6), s = s(r) and ff=ff( 6 0 > 6,<P)- Whereas the other terms are a function of the 
aerosol properties (not aerosol amount) and are contained within the lookup tables.The 
surface reflectance is independent of the aerosol, but dependent on the geometry. In 
practical terms, we parameterize the surface reflectance through the VIS/SWIR surface 
reflectance relationships, which causes it to be a function of scattering angle and 
vegetation index. 

Due to the limited set of aerosol optical properties in the lookup table, the 
equations may not have exact solutions, and solutions may not be unique. Therefore, we 
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find the aerosol solution most closely resembling the set of MODIS measured 
reflectance. In order to reduce the possibility of non-unique retrievals we only allow 
discrete values of rj. Upon completion, the retrieval is assigned a Quality Assurance 
‘confidence’ (QAC) value that ranges from 0 (bad quality) to 3 (good quality). This 
QAC flag is used for creation of Level 3 (gridded) products and for combining land 
retrievals with concurrent over-ocean aerosol retrievals into ‘joint products’ 

5.1 . Selection of “dark pixels” 

Figure 4 illustrates the main steps of the V5.2 land algorithm. Each individual 
MODIS scene, called a granule, consists of a 5-minute swath of data, measuring 
approximately 1340 km by 2030 km. The relevant Level 1 B (LIB) data include 
calibrated spectral reflectance in eight wavelength bands at a variety of spatial 
resolutions, as well as the associated geo-location information. The spectral data include 
the 0.66 and 0.86 pm channels (MODIS channels 1 and 2 at 250 m resolution), the 0.47, 
0.55, 1.24, 1.64 and 2.12 pm channels (channels 3, 4, 5, 6 and 7 at 500 m), and the 1.38 
pm channel (channel 26 at 1 km). The geo-location data are at 1 km and include angles 
( 60 , 6 , 4 ), and 0), latitude, longitude, elevation and date. 0is the scattering angle. The LIB 
reflectance values are corrected for water vapor, ozone, and carbon dioxide (described in 
Remer et al. 2006) before proceeding. 

The first step is to organize the measured reflectance into nominal 10 km by 10 
km boxes (corresponding to 20 by 20, or 40 by 40 pixels, depending on the channel). 

The 400 pixels in the box are evaluated pixel by pixel to identify whether the pixel is 
suitable for aerosol retrieval. Clouds (Martins et al., 2002), snow/ice (Li et al., 2004) and 
inland water bodies (via NDVI tests) are considered not suitable and are discarded. 
Details of this masking are also described in Remer et al. (2006). 

The non-masked pixels are checked for their brightness. Pixels having 2.12 pm 
measured reflectance between 0.01 and 0.25 are grouped and sorted by their visible 
reflectance. The brightest (in the visible) 50% and darkest 20% are discarded, in order to 
reduce cloud and surface contamination and scale towards darker targets. If there are at 
least 12 pixels remaining (10% of 30% of the original 400), then the reflectance in each 
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channel is averaged, yielding the “MODIS -measured” spectral reflectance ff 0A7 , ff 0 . 66 , 
(f 2 .n> and ff l 12 4 - These reflectance values are used for Procedure A. If less then 12 pixels 
remain, then Procedure B (described later) is followed. 

5.2. Correcting the LUT for elevation 

A major change from the old C004 concerns how the algorithm corrects for 
elevated surface targets. The sea-level Rayleigh optical depth (ROD, r RJ ) at a wavelength 
A (in pm) can be approximated over the visible range (e.g. Dutton et al., 1994; Bodhaine 
et al., 1999) by: 

r R<x = 0.00877 A" 405 (12) 

When not at sea level (pressure = 1013 mb), the ROD is a function of pressure (or height, 
z) so that it can be approximated by: 

= Z) = r R Jz = 0)exp(^|) (13) 

where Z is the height (in kilometers) of the surface target and 8.5 km is the exponential 
scale height of the atmosphere. The difference between ROD at z=0 and z=Z is At r x . 

In C004, the algorithm (too) simply corrected the retrieved x product by adding 
the optical depth that was neglected by assuming sea level for the retrieval, (i.e. 
r x {z = Z) = r x (z = 0) + At r x ). However, this correction can give poor results because of 

the large differences between molecular and aerosol phase functions. Instead, the ¥5.2 
algorithm makes use of the procedure described in Fraser et al., (1989). The algorithm 
adjusts the lookup table to simulate different ROD by adjusting the wavelength. 
Substitution of equation 12 into equation 13 yields 

A (z = Z) = A (z = 0)exp( Jj-) . (14) 

For example, at Z = 0.4 km, A increases by about 1.2%. For the blue 0.47 pm 
channel, (centered at 0.466 pm) this means that r R X (z = 0) = 0.194, r R X (z = 0.4) = 0.185 
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and X(z = 0.4) = 0.471 pm. In other words, the algorithm simulates an elevated surface 
by adjusting the blue channel’s wavelength to 0.471 pm. Assuming that gases and 
aerosols are optically well mixed in altitude, the algorithm substitutes for the parameter 
values of the 0.47 pm LUT by interpolating (linearly as functions of log wavelength and 
log parameter) between the 0.47 pm (0.466 pm) and the 0.55 pm (0.553 pm) entries. 
Similar interpolations are performed for the other channels (for example, 0.55 pm would 
be adjusted to 0.559 pm). For the 0.4 km case, this means that lower values of TOA 
atmospheric path reflectance and higher values of transmission are chosen to represent a 
given aerosol model’s optical contribution. However, also note that since the 0.55 pm 
channel has also been adjusted, the associated values of the x indices have been adjusted 
accordingly. 

Whereas most global land surfaces are at sea level or above, a few locations are 
below sea level (Z < 0). In these cases, the algorithm is allowed to extrapolate below 
0.466 pm. Since the extrapolation is at most for a hundred meters or so, this is not 
expected to introduce large errors, and these cases can still be retrieved. Note also that 
due to the extremely low ROD in the 2.12 pm channel, little is gained by adjusting this 
channel. 

5.3. Procedure A: Inversion for dark Surfaces 

If following Procedure A (for dark surfaces), the QAC is initially set to a value 
between 0 (bad quality) and 3 (good quality), depending on the number of dark pixels 
remaining. In Procedure A, the algorithm assigns the fine aerosol model, based on the 
location and time (Levy et al. 2006). From the lookup table, ff, F, T and s (for the fine 
model and coarse model separately) are interpolated for angle, resulting in six values for 
each parameter, each one corresponding to a different aerosol loading (indexed by x at 
0.55 pm). 

The 2.12 pm path reflectance is a non-negligible function of the x, so that the 
surface reflectance is therefore also a function of the x. For discrete values of r] 
between -0.1 and 1.1 (intervals of 0.1), the algorithm attempts to find the x at 0.55 pm 
and the surface reflectance at 2. 12 pm that exactly matches the MODIS measured 


20 



reflectance at 0.47 pm. There will be some error, s, at 0.66 pm. The solution is the one 
where the error at 0.66 pm is minimized. In other words, 

Po.47 “ Po.47 = 0 

Pom - Po.66 = £ (15abc) 

P 2.12 - P 2.12 ~ 0 

where 

P 2 .I 2 ViPz. 12 + FdZ.l 2 P 2 .Ylpl.i 2 10- ~ S {.nPl.u)) + 0 — 7)(p2,12 + Fd £.12^2.^2.12 /(I ~ S 2 .nPz.n)) 

P 0.66 = PiP 0.66 + ^40.66^0.66 f iPl-lz) ^0 — S 0 . 66 f iPz.u))) + 0 “ V)(Po.66 + P ifi. 66 ^ 0 . 66 f iPziz) ^0 “ S 0 . 6 ef iPznff) 

and 

Po.47 = ViPo.4 7 + Fd 0 . 47 T 0 A 78 iP 0 . 66 ) tQ- ~ S a47S (Po.66))) + 0 - h) (Po.47 + 0.47^0.47# (P 0 .& 5 ) — S 0.478(p0.66)))’ 

(16abc) 


where in turn, p a =p a ( r), F=F( r), T=T(r), s = s(r) are functions of t indices in the lookup 
table that is calculated separately for fine and coarse models. Jiff 2 . 12 )’ giffo.6<s ) are 
described by Equations 7-9. Note that non-physical values of rj are tried (-0.1 and 1.1) to 
allow for the possibility of inappropriate assumptions in either aerosol models or surface 
reflectance. Again, the primary products are x 055 , r] a55 , and the surface reflectance (ff zl J). 
The error e is also noted. 

5.4. Procedure B: Alternative Retrieval for Brighter surfaces 

The derivation of aerosol properties is possible when the 2.12 pm reflectance is 
brighter than 0.25, but is expected to be less accurate (Remer et al., 2005), due to 
increasing errors in the VIS/SWIR relationship. However, if Procedure A was not 
possible, but there are at least 12 cloud-screened, non-water pixels, satisfying 

0.25 <p 2 m 12 < 0.25G < 0.40 (17) 

where 
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G-0.5((l//*) + a/V^)) , ( 18 ) 

then Procedure B is attempted. In this relationship pio is cosine of the solar zenith angle, 
cos(0 o ), and ji is cosine of the satellite view angle, cos(0). In procedure B, the QAC is 
automatically set to 0 (“bad quality”). 

Procedure B is analogous to “Path B” described in Remer et al., (2005). Like in 
C004, the Continental aerosol model is assumed. Unlike C004, the VIS/S WIR surface 
reflectance assumptions are those described by Equations 7-9, and the Continental 
aerosol properties are indexed to 0.55 pm. In other words, V5.2 uses equations 10-11, 
except with the first term only (i.e. r\ = 1.0). The primary products for Procedure B are x 
(t 055 ) and the surface reflectance (ff Z i 2 )- The error e is also saved. 


5.5. Derivation of Fine Mode x, Mass Concentration and other 
secondary parameters 


Following the derivation of primary products by Procedure A (t 055 , ? 7 a5J and p s 2 .i^), 
a number of secondary products can also be calculated. These include the fine and coarse 
mode optical depths / ft55 and if ass : 

T 0.55 = ^ 0.5511 0.55 an< ^ 1*0.55 ~ T 0 . 5 s (^ ~ V0.55) ( 1 ^) 

the mass concentration, M: 

+ < 2 «) 

the spectral total and model optical thicknesses t!, and 

t a = t a+< 

where (21) 

r { ~ t o.55^Qx IQ-1. 55 ) an ^ Fa “ t o.55(Q\ (Q 0 . 55 ) 
the Angstrom Exponent a: 

a = ln(r 047 /r 066 )/ln(0.466 /0.644) (22) 
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and the spectral surface reflectance f? K computed by re-arranging Equations 7-9. and 
Af c are mass concentration coefficients for the fine and coarse mode, whereas Q[ x and 
Q\ represent model extinction coefficients at wavelength, A. If the resulting products are 
inconsistent, then the QAC value initially assigned to the pixel is changed to 0 (‘bad 
quality’). If Procedure B was followed, the only secondary products calculated are M and 
t 047 , and the QAC is set to 0. The other products in Procedure B are left undefined. Levy 
et al., (2006) describe how the Q and M c coefficients are defined. 

5.6. Low and negative optical depth retrievals 

A major philosophical change from C004 to C005 is that negative x retrievals are 
allowed. Given that there is both positive and negative noise in the MODIS observations, 
and that surface reflectance and aerosol properties may be under or over-estimated 
depending on the retrieval conditions, it is statistically useful to allow retrieval of 
negative x. In fact it is necessary for creating an unbiased dataset from any instrument. 
Without negative retrievals the x dataset is biased by definition. However, a large 
negative retrieval indicates a situation outside the algorithm’s solution space and should 
not be reported.The trick is to determine the cutoff between a retrieved x that is 
essentially the same as zero, and a retrieved x that is truly wrong. MODIS should retrieve 
with the expected error defined by Equation 1, then values down to -0.05 are essentially 
the same as a zero retrieval and are reported as retrieved. Allowing for slightly higher 
uncertainty, we include x retrievals down to -0. 10 (twice the expected error in pristine 
aerosol conditions), but report these values as -0.05 and lower the QAC value. Note that 
all retrievals With -0.05 < X <0 are reported with high QAC value = 3, unless identified as 
poor quality for some other reason. Some of the products that are retrieved or derived 
(such as r] or Angstrom Exponent) are set to zero or reported as not defined for negative 
retrievals. In cases of low x (x < 0.2), r| is too unstable to be retrieved with any 
accuracy. Therefore, r] is reported as un-defined even though other parameters (such as 
Angstrom exponent and Fine x) may be reported. 
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5.7. Sensitivity study 


Following the lead of Tanre et al ( 1997), we have tested the sensitivity of 
Procedure A by applying it for the following exercises: (1) simulation of conditions that 
are included within the LUT, (2) simulations where one of the parameters (i.e. x) is not 
included within the LUT, and (3) simulations for conditions that include one or more 
errors. 

Exercise 1: Whereas the study of Tanre et al, (1997) tested the algorithm on a 
single geometrical combination, we performed the study in (1) by simulating the 720 
reasonable geometrical combinations in the LUT (0°^^180°, 6$6 0°, 0 O <48°). We 
assumed the “fine” aerosol model to be the neutral (SSA ~ 0.9) aerosol model and that 
the “coarse” model was our Spheroid (dust) model (Levy et al. 2006). For each 
combination of geometry, and for each MODIS channel, we extracted the fine and coarse 
mode values of atmospheric path reflectance p a x , backscattering ratio s x , downward flux 
F d and transmission T x . We assumed that the 2. 1 [xm surface reflectance fS 2 .u = 0- 15, and 
the C004 VIS/SWIR surface reflectance ratios (i.e, ff a6l g= 0-5 p S 2 .i 2 and (fq. 4 f= 0.5 
Using Equations 10- 1 1 , we simulated TOA reflectance p\ for 5 discrete values of r| (r) = 
0.0, 0.25, 0.5, 0.75 and 1.0). Therefore, for each value of x in the LUT, there are 720 x 5 
= 3600 attempts to retrieve that x. 

For smaller x (x ^ 1), the x was retrieved within Ax < 0.01 for all 3600 attempts. 
As x increases, however, computational instabilities lead to a less exact solution. Still, 
though, the retrieved x is certainly within 10% and in most cases to within Ax<0.1. When 
we hold x constant (x = 0.5) and attempt to retrieve values of r] within the LUT On = 0.0 
or 1.0) for the 720 geometrical combinations, we find that rj is retrieved exactly. 

Figs 5 and 6 provide another way of assessing the retrieved MODIS products. Fig. 
5 plots retrieved x, surface reflectance and fitting error as a function of either air mass 
(top) or scattering angle (bottom), given that the input conditions are x OJ!5 =0.5, rp0.5 and 
ffziz= 0.15. In this case, we plotted all of the 720 geometrical combinations in the LUT. 
Air mass is defined as l/cos(9 0 ). The retrieval never exactly matches the input 
reflectances, although the errors are very small (less than 0. 1%). Note that the retrieval 
uses an under-estimated surface reflectance to balance the over-estimated optical depth.. 
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Fortunately, though, most errors are small, and are well within any expected error bars. 
Fig 6 is similar, but for rpO.25, and plotted only for the air mass dependence. The errors 
are much larger (up to 1%), but x is still well within expected error. 

Exercise 2. We used the' same combination of radiative transfer codes (MIEV + 
RT3) used for the LUT (Levy et al., 20G6) to simulate additional values of aerosol 
loading (x 05S = 0.35, 1.5 and 6.0) to create an “extended” LUT. As in exercise (1) we 
simulated the same 720 geometrical combinations as in the V5.2 LUT and the five values 
of T]. On average the retrieval is very close to the expected value, however, the standard 
deviation over all geometry is larger than for x in the normal LUT. A notable exception is 
the attempt at retrieving t 0 55 = 6.0, where the algorithm does a poor job of extrapolating. 
In the operational algorithm, we constrain the maximum possible x to be 5.0. As for 
retrieving values of r\ not included in the V5.2 LUT, the algorithm is successful. The 
r)=0.5 retrieval is well behaved. The attempt at resolving either r)=0.25 or r]=0.75 leads to 
retrieving r|=0.20 and r]=0.70. Although it is impossible for an exact retrieval, due to the 
algorithm choosing between 0.1 intervals, it is interesting that no retrievals of r|=0.30 or 
il= 0.80 are produced. 

Exercise 3. This exercise studied the impact of different types of errors that 
could creep into the retrieval process. Potential errors include (but are not limited to) 
random, systematic or spectrally dependent errors that arise from issues like sensor 
calibration, assuming the wrong aerosol model at a given location, coarse input 
topography mapping, or wrong estimates of the VIS/SWIR surface reflectance 
relationships. These errors are expressed by adding fandom or systematic errors in the 
measurements of one or more spectral channels, geometrical conditions or other input 
boundary conditions. Table 1 lists some prescribed errors, and Table 2 lists eight sample 
geometries used in this exercise. Table 3 shows results when attempting to retrieve 
conditions of t 055 =0.5, r}=0.5 and ff zl2 = 0.15, for the eight sample geometries described 
in Table 2. Table 3A displays the retrieved values of x 0 ^ 5 for each case. Table 3B shows 
the Mean Squared Error (MSB) for each retrieved product, computed over all eight 
geometries. For any case of prescribed errors/geometry, one or more products may be 
over-estimated or under-estimated. If all geometry leads to either one direction or the 
other, the MSE value is designated by (+) or (-). For example, when retrieving with no 
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additional errors (‘LUTinput’), x is never retrieved exactly, but is over-estimated by an 
average MSE of 0.001 1 (+). In balance, ff zl2 is consistently underestimated (MSE of 
0.0004 (-)), with a nonzero fitting error, 8. This is simply a result of computer round off 
error. 

Under most conditions, introducing minor calibration or random errors does not 
destroy the retrieval of x. For most individual errors, the retrieved x is accurate to within 
0.02. However, even when we combine errors (model error, random error, surface error, 
calibration error and geometrical error), we still retrieve x=0.5 with MSE = 0.10, thus 
retrieving within the expected error of Ax=0.125. Retrieval of surface reflectance seems 
to be extremely robust. Retrieval of r] is much more unstable. For simple calibration and 
geometrical errors, the MSE for r) is < 0.1. Combinations of errors lead to large MSE 
(>0.2) for T) retrieval, meaning that ti is not a stable product. Yet, these sensitivity tests 
indicate that generally, the V5.2 aerosol over land algorithm can retrieve useful products. 


6. The aerosol products 

Examples of the three primary aerosol products (x 0 55 , r] and ff zl ^ are shown in 
Fig. 7, along with a color composite of the LIB reflectances (0.47, 0.55 and 0.66 jam 
channels). This image was taken on May 4, 2001 over the U.S. East Coast, and is the 
same image used by King et al., (2003). We note the continuity of the x from land to 
ocean, and that the retrieval of r\ and surface reflectance seem reasonable. Note that rj is 
not plotted over land when x < 0.2. 

Table 4 lists the aerosol over land products that are contained in each “M?D04” 

L2 granule file (MOD04 for Terra and MYD04 for Aqua). For each product, the table 
lists its name within the file, its dimension, and its type. All products are at least two- 
dimensional (nominally 135 x 204 at 10 km x 10 km resolution), and many have three 
dimensions. If there is a third dimension, the channels (usually wavelengths) are listed. A 
parameter’s type may be Retrieved, Derived, Diagnostic, Experimental, or Joint Land and 
Ocean. A Retrieved parameter is one that results directly from the inversion (Procedure 
A), whereas those Derived (such as the Angstrom Exponent), result from those directly 
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retrieved. Products that are Diagnostic include QA parameters and those parameters that 
were calculated during intermediate steps. These diagnostic parameters can be used to 
understand how the retrieval worked. Products denoted Experimental are superfluous to 
the main inversion, may be useful for other applications, but are not discussed here. 
Finally, Joint Land and Ocean products are those that are composites of over-land and 
over-ocean aerosol retrievals. These are intended either for quantitative use (Quality 
Assured where QAC>0; e.g. Optical_Depth_Land_And_Ocean) , or for qualitative 
imaging (QAC s 0; e.g. Image_Optical_Depth_Land_And_Ocean). 

7. Provisional validation of V5.2 products 

The primary means of MODIS validation is by comparing the products with 
equivalent measurements from AERONET or other aerosol measurements. In this way, 
some of the products of C004 (i.e. V4.2 and before) were validated (e.g. Remer et al., 
2005), meaning that they were demonstrated to be accurate to within certain errors. In the 
case of the land products (through V4.2), this meant that ~60 % (slightly less than one 
standard deviation) of the AERONET-measured t values were retrieved by MODIS to 
the expected error described by Equation 1. The other land parameters were either not 
yet validated, or are diagnostic parameters that cannot be validated. 

Since that paper, the algorithm has gone through some minor updates. The last 
update to the C004 family was known as Version 5.1 (‘V5.1’). V5.1 updated the snow 
mask (Li et al., 2005) and cleaned up confusing information in the output files. Even 
though V5.1 never became operational, it is being used in this paper to bridge between 
the operational C004 family of algorithms and the new algorithm we introduce in this 
paper. In addition to validation by AERONET, we make qualitative analyses based on 
visual inspection and global statistics. 

7.1 . Direct comparison of V5.2 and V5.1 products 

Fig. 8 plots retrieved x at 0.55 pm from both V 5.1 and V5.2, over small areas of a 
MODIS granule. V5.1 (OLD) is presented in (a), whereas V5.2 (NEW) is shown in (b). 
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Fig. 8 shows a region in the western U.S. from 30 Sep 2003. The V5.2 aerosol retrieval 
adds more valid retrievals over very low x areas (coastal Oregon and northern California). 
V5.2 reports these areas as having near zero or slightly negative x, where V5.1 would 
have reported fill values (errors). In areas farther from the coastline, ¥5.2 tends to clean 
up contamination presumably caused by clouds, elevation, and inhomogeneous surface 
properties, and produces a much more reasonable picture of x. 

7.2. Statistics of V5.2 versus V5.1 

Of most interest to the climate community will be the changes in the statistics of 
the aerosol products. These include the global mean values and the distribution 
(histogram) of the values. For the set of MODIS granules listed in Table 5 (about 6300 
granules of both Terra and Aqua), the mean 0.55 pm x is reduced from 0.28 to 0.21. This 
is a significant reduction that can be compared with model estimates. 

Fig 9 plots the histograms of retrieved x at 0,55 pm from both Y5.1 and V5.2. 
These histograms include 141 individual Terra and Aqua granules that are known as the 
MODIS “test_bed”, and twelve days of global data - all listed in Table 5. The use of 
global data is especially important for determining how the retrieval behaves in regions 
not selected for algorithm development. Of course, the obvious change in the V5.2 
product is that small magnitude negative x retrievals are valid. About 10-11% of the total 
x retrievals are now retrieved as below zero, of which only about 3% are below -0.05. 

This promising result indicates that V5.2 has reasonable ability to detect very clean 
conditions within the expected error of ±0.05. Also noted in Fig 9 is that the fraction of 
retrieved medium to medium high x (0.2 < x < 0.75) is reduced, while the fraction of high 
x (x > 0.75) remains about the same. 


7.3. Comparison of V5.2 to V5.1 and with AERONET 



As of 1 April 2006, the V5.2 algorithm has been run on nearly 6300 granules, 
including one full month (August 2001), fifteen entire days (listed in Table 5) and about 
141 individual granules that are known as the MODIS “test_bed.” These granules include 
observations from both Terra and Aqua, and are seasonally and yearly representative of 
the MODIS time series. For comparison, we ran V5.1 on the same set of granules. 

Figures 10 and 1 1 plot the comparisons of both V5.1 and V5.2 with the AERONET data, 
via the spatio-temporal co-location method of Ichoku et al., (2002). 

Figure 10 plots the retrieved MODIS x against AERONET x, both at 0.55 pm. 

The data have been sorted by AERONET x and averaged into bins with equal numbers of 
observations in each bin. The mean and standard deviation of each bin are calculated and 
plotted in Fig. 12 as a solid dot and error bars. The correlation is calculated from the 
freely plotted points before binning, although the cloud of points is not shown in the plot. 
The regression equation has improved tremendously, from “y=0.097+0.91x” to 
“y=0.029+1.01x.” Correlation R is also improved, from R=0.847 to R=0.894. It should 
be noted that slight differences in the number of points arise due to different selection of 
valid dark pixels and allowance of below zero x retrievals. 

Figure 11a plots MODIS q against AERONET rj, where AERONET q is 
calculated from sun observations of spectral x as described by O’Neill et al., (2003). 

Keep in mind that unlike MODI S/AERONET comparisons of x, MODIS and AERONET 
do not retrieve the same quantity labeled as r|. The AERONET retrieval assumes one 
fine mode and one coarse mode. Thus, AERONET r| is the weighting between modes. 
The MODIS land q is a weighting between models, where fine-dominated models 
contained coarse modes and coarse-dominated models contain fine modes. The 
improvement to the MODIS q product is mainly its correlation to AERONET. Note that 
q is defined only when x > 0.2. Figs lib and e show comparisons for derived products, 
including the Angstrom Exponent (defined by 0.47 and 0.66 pm), and Fine optical depth 
(i.e. x f = x x q), respectively. For Fine x, the correlation and slopes are nearly unchanged 
between V5. 1 and V5.2; however, the offset goes from +0.05 1 to -0.031. The result is 
that nearly two-thirds of all V5.2 MODIS Fine x fall within expected errors defined by 
Eq. 1. Note again that the difference in the number of points is due to different selection 
of dark pixels and treatment of negative x retrievals. The Angstrom exponent has little 
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improvement from V5. 1 to V5.2, except for slightly better but still poor correlation with 
the AERONET measured quantities. In general, the changes to the MODIS aerosol 
retrieval algorithm described here have resulted in a much less biased x and x f products 
than the previous algorithm. MODIS r\ correlates better with AERONET, although it still 
leaves room for improvement. 

8. Conclusion 

In this document, we have introduced a new algorithm (V5.2) for deriving aerosol 
optical properties over dark land surfaces, from MODIS observed spectral reflectance. In 
the new algorithm, we have updated a number of assumptions, including the VIS/S WIR 
surface reflectance parameterization, and the statistical implications of deriving below 
zero aerosol optical thickness. Most significantly, instead of an independent two-channel 
retrieval, V5.2 is a simultaneous three-channel inversion that makes use of aerosol 
information contained in the SWIR (2.12 pm) channel. We have coupled these changes 
with updated representative global aerosol optical models that are described in Levy et al. 
(2006). 

The V5.2 algorithm has been tested, both for its theoretical ability to derive 
aerosol properties, and on a test bed of 6300 MODIS granules. Compared with co-located 
AERONET sites, the V5.2 MODIS algorithm retrieves aerosol properties more accurately 
than V5.1. Specifically, the retrievals of total x meet expected accuracy levels 
(±0.05±0.15x). MODIS/ AERONET x regression has an equation of: y = l.Olx + 0.03, R 
= 0.90. Global (the 6300 granules) mean x has been reduced from 0.28 to 0.21. 

Retrievals of r) show less significant improvement, but are still better correlated with 
AERONET results than previous versions. Retrievals of spectral Angstrom Exponent 
show little or no improvement at this time. However, the new algorithm’s derivation of 
Fine x (x x r|) is much improved. This product that can be related to the anthropogenic 
contribution to the total x (e.g. Kaufman et al., 2005) and has specific applications for the 
air quality community. Finally, the V5.2 products’ quality assurance (QA) has been 
overhauled and is now more useful to users within the aerosol community. 
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TABLE 2: Solar/Stjrface/Satellite Geometry for eight examples 


Reference 

Solar Zenith 

View Zenith 

Relative Azimuth 

Scattering Angle 

A 

12.00 

6.97 

60.00 

163.40 

B 

12.00 

52,84 

60.00 

120.53 

G 

12.00 

6.97 

120.00 

169.59 

D 

12.00 

52.84 

120.00 

132.35 

E 

36.00 

6.97 

60.00 

140.12 

F 

36.00 

52.84 

60.00 

104.74 

G 

36.00 

6.97 

120.00 

147.00 

H 

36.00 

52.84 

120.00 

136.29 


All units are degrees 
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TABLE 1: LIST OF PRESCRIBED ERRORS FOR V5.2 SENSITIVITY STUDY 


Reference 

Error Name 

Description 

1 

LUTinput 

LUT input: Use the LUT with no prescribed errors 

2 

ModError 

Aerosol model error: We tried to retrieve with the Non-absorbing fine model LUT 

3 

RndError 

Random Error: All channels have random reflectance error of up to ±0.002 

4 

SfcError 

Surface Error: 10% error in assumed 0.66/2.12 surface reflectance relationship 

5 

CalError 

Calibration Error: All channels have random error of up to ±1% 

6 

ElvError 

Elevation Error: Elevation is 1km instead of assumed sea level 

7 

GeoError 

Geometry Error: All angles have random error of up to ±5 degrees 

8 

AllError 

Combination of 2,3 ,4,5,6 and 7. 
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TABLE 3: Results of Sensitivity Study using Prescribed Errors 


Geometry 
Error Name 

LUTiupiit 

RadError 

CalError 

GeoError 

ModEiror 

ElvError 

SfcError 

AIlError 

A 

0.501 

0.4786 

0.5242 

0.5143 

0.5015 

0.6068 

0.5402 

0.6963 

B 

0,501 

0.4887 

0.5242 

0.4977 

0.4993 

0.6035 

0.5422 

0.6677 

C 

0.501 

0.5227 

0.5227 

0.4657 

0.4835 

0.5104 

0.4955 

0.4809 

D 

0.5011 

0.5104 

0.4995 

0.4761 

0.5014 

0.5228 

0.498 

0.4892 

E 

0.5008 

0.4754 

0.502 

0.4893 

0.4866 

0.5211 

0.4877 

0.5737 

F 

0.501 

0.5135 

0.5029 

0.4922 

0.5035 

0.531 

0.488 

0.5536 

G 

0.5014 

0.4973 

0.5199 

0.4698 

0.4811 

0.5097 

0.488 

0.427 

H 

0.5016 

0.4961 

0.5001 

0.4744 

0.5198 

0.5299 

0.4939 

0.5106 


A: Retrieved x at 0.55 \im (expected x=0.5) 


Product T T7rri . , 

Error Name LUTmput 

RndError 

CalError 

GeoError 

ModError 

ElvError 

SfcError 

AIlError 

T 

0.001 1(+) 

0.0159 

0.0162 

0.0215 

0.0123 

0.056 1(+) 

0.0221 

0.1006 


0.0000 

0.0000 

0.0707 

0.1000 

0.0707 

0.4243 (+) 

0.1323 (+) 

0.4912 (+) 

p 

0.0004 (-) 

0.0008 

0.0022 

0.0025 

0,0031 (-) 

0.0067 

0.0020 (+) 

0,0074 (+) 

E 

0.0010 

0.0021 

0.0037 

0.0028 

0.0020 

0.0025 

0.0035 

0.0052 


B: MSE of retrieved x,r),p° and e (expected t= 0.5, T)=Q.5, p°=0.15 and e=0.0). Entries designated with (+) mean that the product 
was over-estimated for all 8 geometries, whereas those with a (-) means it was under-estimated for all geometries. 
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TABLE 4: Contents of MODIS Y5.2 Aerosol Level 2 File (MQD04/MYD04): Land Products 


Name of Product (SDS) 

Dimesions : 3 rd Dimension 

Type of product 

Corrected_Optical_DepthJLand 

X,Y,3: 0.47, 0.55, 0.66 ptm 

Retrieved Primary 

Corrected_Optical_Depth_Land_wav2p 1 

X,Y,1: 2.12 \im 

Retrieved Primary 

Optical„DepthJRatio_SmallJLand 

X,Y: (for 0.5:5. pm) 

Retrieved Primary 

Surface_Reflectance_Land 

X,Y,3: 0.47, 0.66,2.12 \m 

Retrieved Primary 

FittingJError_Land 

X,Y: (at 0.66 jxm) 

Retrieved By-Product 

Quality_Assurance„Land 

X,Y,5: 5 bytes 

Diagnostic 

Aerosol_T ype_Land 

X,Y: 

Diagnostic 

Angstrom„Exponent_Land 

X,Y: (for 0.66/0.47 \xm) 

Derived 

Mass_Concentration„Land 

X,Y: 

Derived 

Optical_Depth_Small_Land 

X,Y,4: 0.47,0.55,0.66,2.12 pm 

Derived 

Mean_Reflectance_Land 

X,Y,7: 0.47, 0.55, 0.66,0.86,1.2,1.6, 2.12|xm 

Diagnostic 

STD_Reflectance_Land 

X,Y,7: 0.47, 0.55 ,0.66,0.86, 1.2,1. 6,2. 12|xm 

Diagnostic 

Cloud_Fraction_Land 

X,Y: 

Diagnostic 

Number_Pixels_Used_Land 

X,Y: 

Diagnostic 

Path_Radiance_Land 

X,Y,2: 0.47, 0.66 |xm 

Experimental 

Error_Path_Radiance„Land 

X,Y,2: 0.47, 0.66 fxm 

Experimental 

Critical_Reflectance_Land 

X,Y,2: 0.47, 0.66 pm 

Experimental 

Error_Crit_Reflectance_land 

X,Y,2: 0.47, 0.66 pm 

Experimental 

Error Critical JReflectance_Land 

X,Y,2: 0.47, 0.66 pm 

Experimental 

Quality _Weight_Path_Radiance_Land 

X,Y,2: 0.47, 0.66 pm 

Experimental 

Quality_Weight_Crit„Reflectance_Land 

X,Y,2: 0.47, 0.66 pm 

Experimental 

Optical_Depth_Land_And_Ocean 

X,Y: 0.55 pm 

Joint Land and Ocean 

hnageJ3pticalJDepth_Land_And_Ocean 

X,Y: 0.55 pm 

Joint Land and Ocean 

Optical_Depth_Ratio_Small_Land_And_Ocean 

X,Y: 0.55 pm 

Joint Land and Ocean 


X = 135; Y = 203. If there is a 3 rd dimension of the SDS, then the indices of it are given. The "Retrieved” parameters are the 
solution to the inversion, whereas “Derived” parameters follow from the choice of solution. “Diagnostic” parameters aid in 
understanding of the directly Retrieved or Derived products. “Experimental” products are unrelated to the inversion but may have 
future applications. “Joint Land and Ocean” indicate combined land and ocean products. 
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Atmospherically Corrected Surface Reflectance 
Visible versus 2.1 pm 




Suifaee Reflectance <0.644 i*m) 

Figure 1: Atmospherically corrected surface reflectance in the visible (0.47 and 0.66 
pm channels) compared with that in the 2.12 pm SWIR channel (a), and the 0.47pm 
compared with that in the 0.66 pm channel (b). 
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Surface Reflectance Relationship vs Scattering Angle 
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Surface Reftactance Relationship vs Scattering Angle 
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Figure 2: VIS/S WIR surface reflectance relationships as a function of scattering 
angle. The data were sorted according to scattering angle and put into 20 groups of 
equal size (about 230 points for each scattering angle bin). On all subplots, each 
point is plotted for the median value of scattering angle in the bin. Part (a) plots 
median values of reflectance at each channel as a function of the scattering angle. 
Linear regression was calculated for the 230 points in each group. The slope of the 
regression (for each angle bin) is plotted in (b), the y-intercept is plotted in (c) and 
the regression correlation is plotted in (d). Note for (b), (c) and (d) that 0.47 pm vs 
2.12 pm (r0470) is plotted in blue, 0.66 pm vs 2.12 pm (r0660) is plotted in red and 
0.47 vs 0.66 pm (rvis) is plotted in green. 
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0.644 p.m VS 2.12 j nn Surface Reflectance 
as a function of NDVI bins 

SWIR 



Surface Reflectance (2.12 pm) 

Figure 3: 0.66 pm versus 2.12 pm surface reflectance as a function of bins of 
NDVI swm values. The standard regression is plotted, with regression equations given 
in the lower right hand corner. The ratios (if forced through zero) are given 
beneath the legend. Blue refers to low NDVl^g, red to medium and green to high 
values. 


.43 




MOOIS Over Land Algorithm 

All procedures applied to individual boxes of 20x 20 pixels at 500 m resolution (10 km at nadir) 
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Figure 4: Flowchart illustrating the derivation of aerosol over land for V5.2 
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Retrieved x vs Scattering Angie 
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Figure 5: Retrieved MODIS products as a function of Air Mass (a-c) and Scattering 
Angle (d-f) for inputted atmospheric conditions (t=0.5, ri=0.5 and p s 2 12 =0.15) and 
720 LUT geometrical combinations. The retrieved x is plotted in (a) and (d), the 2.12 
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pm surface reflectance in (b) and (e) and the fitting error is plotted in (c) and (f). 
Note that in all cases, the rj value of 0.5 was retrieved exactly. 
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Figure 6: Retrieved MODIS products as a function of Air Mass for inputted 
atmospheric conditions (t= 0.5, r|=0.25 and p S 2. 12 =0.15) and 720 LUT geometrical 
combinations. The retrieved r is plotted in (a), retrieved T| in (c), the 2.12 pm 
surface reflectance in (c) and the fitting error is plotted in (d). 
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A) RGB: p M {0.66, 0.55, 0.47 pm) B) -t(0.S5 pm) 



C) r|(0.55 pm) D) p s (2.1 pm) 


Figure 7: Retrieved aerosol and surface properties over the Eastern U.S. on May 4, 
2001. This figure can be compared with that plotted in King et al., (2003). Panel A) 
is a ‘true-color’ composite image of three visible channels, showing haze over the 
mid- Atlantic. Panels B) and C) show retrieved t and rj, showing that the heavy 
aerosol (t ~ 1.0) is dominated by fine particles. The transport of the aerosol into the 
Atlantic is well represented with good agreement between land and ocean. In fact 
the continuity of t seems to be improved since earlier versions of the aerosol 
algorithm. Note that over-land rj is not reported when x < 0.2. Panel D) shows the 
retrieved surface reflectance. 
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Sep 30, 2003; 17:55 UTC 



Figure 8: Retrieved t (AOT) at 0.55 p,m for Old V5.1 (a) and New V5.2 (b) over 
California for 30 September 2003. The color scale is the same for both plots. Note 
the increase in the retrieval spatial coverage and reduction in surface contamination 


for V5.2. 
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Figure 9: Histogram of retrieved t (AOD) over land, from V5.2 (C005) in green, 
compared to V5.1 (C004) in orange. The data include the 141 granules of the Terra 
and Aqua “testjbed” as well as twelve complete days. The value of each bin refers to 
the minimum value of the bin (the max value would be the value of the next bin). 
Note that the general lognormal nature of the retrievals is preserved, except now 
there are some negative values. 
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AERONET x at 0*g5pm {Quadratic Fit) 

Figure 10: MODIS t over land retrieved at 0.55 pm, compared with AERONET t 
interpolated to 0.55 pm. The solid shapes and error bars represent the mean and 
standard deviation of the MODIS retrievals, in 20 bins of AERONET -derived t. 
Both the retrievals from V5.1 (orange) and V5.2 (green) are shown. The regressions 
(solid lines) are for the cloud of all points before binning (not shown). The expected 
errors for MODIS (±0.05 ±0.15t) are also shown (dashed lines). 
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Figure 11: MODIS aerosol size retrievals compared with AERONET derived 
products. The solid shapes and error bars represent the mean and standard 
deviation of the MODIS retrievals, in 20 bins of AERONET -derived product Both 
the retrievals from Y5.1 (orange) and V5.2 (green) are shown. The regressions (solid 
lines) are for the cloud of all points (not shown). A) rj over land retrieved at 0.55 
pm, compared with AERONET t] retrieved by the O’Neill method. Note that r| is 
defined differently for MODIS and AERONET and that we only show results for x 
> 0.20. B) MODIS-derived a (0.466/0.644 pm) over land with AERONET a 
interpolated to the same wavelengths. O MODIS Fine x over land retrieved at 0.55 
pm, compared with AERONET Fine x interpolated to 0.55 pm by quadratic fitting 
and the O’Neill method. The expected errors for MODIS (±0.05 ±0.15x) are also 
shown (dashed lines). 
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